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We show how classical and quantum dualities, as well as duality relations that appear only in a 
sector of certain theories {emergent dualities), can be unveiled, and systematically established. Our 
method relies on the use of morphisms of the bond algebra of a quantum Hamiltonian. Dualities 
are characterized as unitary mappings implementing such morphisms, whose even powers become 
symmetries of the quantum problem. Dual variables -which were guessed in the past- can be derived 
in our formalism. We obtain new self-dualities for four-dimensional Abelian gauge field theories. 

PACS numbers: 03.65.Fd, 05.50.+q, 05.30.-d 



Introduction. Dualities appear in nearly all disciplines 
of physics and play a central role in statistical mechanics 
and field theory [l|, Q. When available, these mathe- 
matical transformations provide an elegant, efficient way 
to obtain information about models that need not be 
exactly solvable. Most notably, dualities may be used 
to determine features of phase diagrams such as bound- 
aries between phases, and the exact location of some 
critical/multicritical points. Historically, dualities were 
introduced in classical statistical mechanics by Kramers 
and Wannier (KW) as a relation between the partition 
function of one system at high temperature (or weak cou- 
pling) to the partition function of another (dual) system 
at low temperatures (or strong coupling). This relation 
allowed for a determination of the exact critical temper- 
ature of the two-dimensional Ising model on a square 
lattice 0], before the exact solution of the model was 
available. Later on, it was noticed that, due to the con- 
nection between quantum theories in d space dimensions 
and classical statistical systems in d + 1 dimensions, du- 
alities can provide relations between quantum theories in 
the strong coupling and weak coupling regimes [if. The 
current work is motivated by a quest for a simple unifying 
framework for the detection and treatment of dualities. 

We will describe an algebraic approach to dualities and 
self-dualities for systems of arbitrary spatial dimensional- 
ity d. We will show that quantum (sclf-)dualities (a con- 
nection between Hamiltonians) become dualities of the 
related classical statistical problem in d + 1 dimensions. 
Thus, quantum and classical (self-) dualities are intrin- 
sically equivalent, yet it will become clear that quan- 
tum (self-) dualities are -with the technique presented 
here- much easier to detect and exploit. The gist of the 
method is the characterization of quantum (self-)dualitics 
as structure preserving mappings (homomorphisms) be- 
tween operator algebras which are Hamiltonian depen- 
dent. The structure of quantum mechanics further re- 
quires that these (self-)duality mappings should be uni- 
tarily implementable. In contrast, generalized Jordan- 
Wigner transformations [1] for example, are dictionaries 
connecting representations, independent of the structure 



of any particular Hamiltonian. 

Bond Algebras and Dualities. Our main thesis is that 
quantum dualities (self-dualities) are homomorphisms 
(automorphisms) of bond algebras [H that preserve lo- 
cality of interactions and can be implemented through a 
unitary map. Take a quantum Hamiltonian H, given as 
a sum of quasi-local operators or bonds {h^\ weighed by 
couplings an, H — ^2 R anIiR. The index R can repre- 
sent, for example, lattice sites. The bond algebra of H, 
Ah, is the smallest operator algebra that contains every 
bond in H -and thus H itself. It can be described as the 
algebra of all linear combinations of products of bonds 
[\ Iir and the identity operator. The core idea is that two 
Hamiltonians H and H^ ua \ are dual to each other if there 
is a unitarily implementable homomorphism $ between 
their bond algebras mapping H to -ffduai up to irrelevant 
terms in the thermodynamic limit [|| . So we demand that 
$(iJ") = UdHUJj = ifduai + Vb where the boundary op- 
erator Vb is irrelevant [6( . If H and -f/duai share the same 
bonds but with different couplings, then the duality is 
nothing but a self-duality, established through an auto- 
morphism of Ah- This scenario includes the very useful 
special case of two exchanged couplings representing a 
weak coupling^strong coupling exchange. To make 
clear that this approach is physically sensible, it is enough 
to notice that such homomorphisms preserve the Heisen- 
berg equations of motion. Notice that the labels {R} are 
completely arbitrary, no reference is made to any partic- 
ular geometry or dimensionality. The primary algebraic 
objects are the bonds (Hj, built out of elementary degrees 
of freedom such as spins. In the past, quantum dualities 
such as KW were presented as non-local mappings be- 
tween elementary degrees of freedom. In contrast, duality 
morphisms arc mappings local in the bonds and, remark- 
ably, provide means to derive those non-local mappings 
(which shows that these self-duality automorphisms are 
indeed the quantum version of the classical order-disorder 
transformations of Kadanoff and Ceva Q ) ■ That all du- 
alities are manifestations of bond algebraic morphisms is 
not obvious. If, however, as is the standard case, two 
systems are dual to one another on general subsets A 
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of an infinite lattice then an exact duality between the 
two systems exists if and only if the bond algebras of 
the two systems are identical. The proof of this asser- 
tion is straightforward. The proviso of general sublattices 
implies that a unitary transformation giving rise to the 
same spectrum may be applied for a general collection 
of bonds R G A and their duals R' G A': U D HU D = 

u o( K J2ReA a R h R) U D = J2n'eA' a R' h 'w = H ^\- As 
this holds for all A, it follows that U ahull D = h' R , for 
all R, R' . If two sets of operators (including the bond 
operators {Iir}) are related by a unitary transformation 
Ud then their algebras are identical. Similarly, if two sets 
of operators {/i_r} and {h' R ,} exhibit an identical algebra 
then there is a unitary transformation Ud relating them. 

In general, self-dualities do not leave H invariant. They 
are symmetries of the bond algebra Ah, and this is the 
key to detect them. However, they may become symme- 
tries on appropriate regions of parameter space. If, e.g., 
Ud exchanges the couplings g and g' in H then at the 
self-dual point g = g' , [H,Ud] = (up to the irrelevant 
terms @). Moreover, if Ud effects the exchange for any 
values of g and <?', then for even n, [U^,H] = (again, 
up to irrelevant terms). Taking n = 2 we see that 



Self-duality — > ^/Quantum Symmetry . 



Thus a self-duality could reveal nontrivial hidden sym- 
metries of a problem. Of course, the symmetries U^, 
need not be all independent or non-trivial (we will see 
examples below). One can always add to H an irrelevant 
boundary term Vg (related, but not equal to Vb) derived 
from the bond algebra, so that even for finite systems 
\H + Vg, Up] = exactly. Thus, it may be useful to work 
with the more symmetric H + Vg . 

As a basic illustration, take H[j, h] = j X^i 1 °f °f+i + 

hY,i=i vf+jvN = H U> h ]+jo z N , (jo z N = V&) (the of are 
Pauli matrices), where H[j,h] is the Hamiltonian of an 
Ising chain in a transverse magnetic field (TV spins). One 
can check that of i-> o z N , o z N i-> of , of i-> o z ^o z ^ i)+1 
(i = 2, 3, • • • , N), o z o z +1 1 ^ <t? w (i = 1, • • • , N-l), with 
r(i) = N+l—i, gives a unitarily implcmentable automor- 
phism $ of H's bond algebra. <!> is clearly a self-duality 
for the Ising chain H[j, h], U D H\j, h]U D = H[h,j] + V B , 
with boundary term Vb = ho z N — jof, and it is an 
exact self-duality for H, &(H\j,h]) = U D H[j,h]U D = 
H[h,j]. In this simple case, U D = 1. The standard 
approach Q to this self-duality involves defining non- 
local spin operators -the dual variables- but nothing in 
principle determines their form; dual variables have to 
be guessed. In contrast, in our formalism it is natu- 
ral to use the duality mapping to define dual variables 
fif as /if = UDofU D . Then the above relations lead 
to /Ltf = o z N , nf = K(i) a r(i)+V i = ,N. On the 

other hand, fxf = U D o z U D = U D o z o* +1 x • ■ • o%_ x o% x 



o z N U D , so that, by the duality mapping above, reduces to 
M? = IlL^V) = n!l! + '^' Similarly, the Jordan- 
Wigner dictionary Q gives rise to a bond algebra map- 
ping when applied to d = 1 spin and spinless Fermi 
systems. The explicit exchange statistics transformation 
can be derived by solving for one set of bonds in terms 
of the other. It can be shown that there is no Jordan- 
Wigncr transformation that relates two local Hamiltoni- 
ans in dimensions d > 1: By examining the product of 
bonds around closed loops an inconsistency is found if 
local spin-less Fermi bilincars could be mapped to local 
spin terms and vice versa. In the following we disregard 
boundary terms without further comments. 

Dualities and Self- dualities in Quantum Statistical Me- 
chanics. The d = 3 orbital compass (OC) model 

^OC = - '^2\JxSjSj +gi + J y S^S^ +S2 + JzSjSi + g 3 ] 

T 

(Sg = \oj) has been proposed [§] to study orbital or- 
dering in transition metal compounds. A still interesting 
yet simplified scenario for orbital ordering is provided by 
the planar OC model (POC) 

Jy] = - E^^+el + J V°M + e 2 ) (1) 
1 

Its bond algebra Ah PO c i s generated by 
{ <J r a T+s 1 ' a T a T+e 2 }^ and it: is specified by a few 
relations: Each bond (i) squares to one, (ii) anti- 
commutcs with the four other bonds which share 
any of its vertices, and (iii) commutes with all other 
bonds. The mapping <P{o§of +Si ) = ol^ol^^ 

H^^7+S 2 ) = ^r+e^r+e 2+e v P^serves every relation 
among bonds, showing a self-duality under J T «-> J y . 
The POC Hamiltonian is dual as well [l0( to the 
Xu-Moore (XM) Hamiltonian [n| 

#xm [j, h] = -J2(jno*?+ ho x r ), (2) 

(with Do z T = o§ol +gi o§ +gi+g2 oS +g2 ) which was intro- 
duced as a simplified model for some aspects of quantum 
phase transitions in p + ip superconducting arrays. The 
duality comes from the mapping of bonds &{&?&? + g 1 ) = 
□cr 2 ,-, $(ofof + g ) = of + g 2 , which is indeed given by a uni- 
tary U D - Thus U D Hpoc[J x , J y ]U D = Hxm[J x , Jy), and 
these two models must have the same phase diagram. In 
spite of this, the quantum(d)-to-classical((f-|- 1) mapping 
is much easier for Hxm than for -ffpoc, another mani- 
festation of the power of duality transformations and a 
useful fact if one wants to perform, say, quantum Monte 
Carlo simulations. The self-duality of the XM Hamilto- 
nian can be deduced from the self-duality of the POC 
model and the duality just described, or directly as an 
automorphism of its bond algebra. Applied to the ele- 
mentary degrees of freedom {cr5, cri}, the automorphism 
returns the non-local dual operators of [ill ]. 
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Classical from Quantum Dualities. The standard 
quantum((i)-to-classical((i + 1) connection establishes 
an equivalence between quantum (as unitary map- 
pings) and classical dualities. Take for example 
the XM Hamiltonian H X u[j,h] of Eq. © 

classical rendition is (| sinh(2J*)) !1 ^ 2 Z[J, K] 

Z[J,K] = Y j[a] e T -^ jn,J ^ {t)+Ka ^ t+1 \ J = 

J* = h^, and K = -± lntanh (/if^ 

length along the time axis N t ;§> 1. 

, n/2 
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•At 

J N t ' 
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larly, H XM [h,j] maps to (± sinh(2J))"^ Z[J* , K 
K* = — i lntanh (jj^j ■ ^ follows already that 
sinh2Jsinh2A'* = 1 = sinh2J* sinh2AT, yet nothing 
in principle guarantees any relation between Z [J, K] 
and Z[J*,K*] so far. Now, due to the quantum 
self-duality HxN\[j,h] = UdHxm[1i, j]Up, we have that 
Trexp (— AtHxm [j, h]) = Trexp (— ArHxM[h, j]). Hence 

which is indeed the 



Z[J,K] 



Z[J*,K*] 



n/2 j 



[1 11 ] by considerably 



(±sinh(2J)J ( ± sinh(2J*)J 

classical self-duality obtained in 
more laborious classical methods. 

Emergent Dualities. A (sclf-)duality can emerge in a 
sector of a theory (e.g., for particular subsets of cou- 
plings, or low energy subspace). The projection of a bond 
algebra onto a sector of the full Hilbert space generates 
a new bond algebra that may have (self-)dualities not 
present in the full model. An example is provided by the 
Quantum Dimcr Model (QDM) [12j defined on the or- 
thonormal set of dense dimcr coverings of a lattice. The 
QDM Hamiltonian reads 



h QD m = EHOnx: 

□ 

+ «(li i) <i il + 



:><nl 



:><~D] 



(3) 



with the sum performed over all elementary plaquettcs. 
The QDM contains both a kinetic (t) term that flips one 
dimcr tiling of any plaquctte to another (a horizontal cov- 
ering to a vertical one and vice versa), and a potential (i>) 
term. At the (so-called) RK point t = v [12], the ground 
states are equal amplitude superpositions of dimcr cov- 
erings. If P g is the projection operator onto the ground 
state sector, then P fl [(|J J) (~ | + |~) (U \)h p a = 



^[(IIIXIII 



!)(nl)b p s = x n p 9' with 



xu = or 1 on the particular plaquette □ where 

[(II IXHI + 1~) (I !!)]□ fli P s thc dimcr in thc p la - 

quette □. At the RK point, the projected Hamiltonian 
becomes P g PtQDMP g = 0. Since both the kinetic (t) 
and potential (v) terms are given by x\jP g within the 
ground state sector, the kinetic and potential operators 
can be interchanged without affecting the bond algebra. 
This self-duality emerges exclusively in the ground state 
sector of the QDM at the RK point. 

Dualities in Quantum Field Theory (QFT). An ele- 
mentary application of our technique is provided by a 



free massless scalar field in 1 + 1 dimensions with 
Hamiltonian H = i J dx [ir 2 (x,t) + ^^gj— ^ ], and 
[4>(x, t), tt(x' , t)] = iS(x — x'). (With obvious modifica- 
tions, this Hamiltonian describes a taut string.) To study 
this model's bond algebra, it is convenient to discretize 
it, with lattice spacing a, i.e. | EJ 71 "* + ~ 4>i) 2 1 ?}- 
The automorphism 7Tj i— > — (<A;+i — <f>i)ja, (<A;+i — 
<j>i)/a i — ► — 7Ti + i preserves the canonical commutation re- 
lations. Thc dual variables provide a convenient way 
to study this self-duality in the continuum. Their dis- 
crete form is < 



m—i+l 



7Tj = -(4> t +i - 4>i)/a. 
Now we can let a go to zero to obtain dual variables 
in thc continuum: Tc(x,t) = Uo'n'{x,t)Ul ) = —^(x,t), 

4>(x,t) = U D <l)(x,t)W D = dy ir(y,t). These are 
toy examples of solitonic variables. In general, self- 
dualities can be destroyed by coupling the system to 
sources, but this is not necessarily the case. Consider the 
scalar field now coupled to external classical sources A, E: 

H A - E = Jdx 



±(.-AA)2 + !(§§) - XE4> 



The self- 



duality maps H ' to 



H 



A,E_ 



= Idx 



-(it - XA) 2 + - ( ^) - \E(f> + c- number 
2 2\axJ 



The self-duality survives this coupling to ex- 
ternal sources, with dual sources A(x,t) = 
J^dw E(w,t), E(x,t) = ^(x,t). 

Next we consider 1m gauge field theories (GFTs) de- 
fined on a Euclidean 3 + 1-dimcnsional lattice. The inter- 
est in these theories grew out of 't Hooft studies on quark 
(charge) confinement in pure SU(N) gauge theories fl3j |. 
that suggest that their most important degrees of free- 
dom near a confincmcnt-deconfinement phase transition 
are the field configurations taking values in the center 
subgroup of SU (A), Zjv- To explore this scenario, several 
author considered Wilson's action for Euclidean lattice 
GFTs Q, S = -£(En Rc TriUijU^UkjU^ - 1)), re- 
stricting the fields to take values in Zjy. This is the model 
we are going to study, thus Ui.j stands for a Ath root of 
unity attached to the oriented link with endpoints i,j, 
and Uij = U* i . In the axial gauge the action simplifies 



S 



1 

25 



2 EE 

i=l 



'n+e 4 



J + cos(e*j] , 



(6 r \ = 9 3 n+e2 -dl- e 2 n+e3 + el), and cyclic permuta- 
tions thereoff. The goal is to learn about duality proper- 
ties of amplitudes in QFTs, as given by a path integral 
over field configurations. Computation of a vacuum to 
vacuum amplitude (0out|0in) amounts to evaluating a 
partition function. Thus we can apply the bond alge- 
bra technique to look for self-dualities in QFTs that are 
more conveniently quantized through path integrals. To 
proceed, we need to compute thc quantum Hamiltonian 
equivalent to thc gauge fixed action given above. This 
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is a difficult task for arbitrary N, but the computations 
were done (in a different context) in 11511 . Using these 
(the coupling K depends on TV and g |16j ) 



A91 



h.c. 



where A9 3 = U^U 2 +ei U n \ e2 U 2 ^ , and cyclic permuta- 
tions. There are now N x N unitary matrices U, V on 
each link (n, ej), i = 1, 2, 3 of a cubic lattice, (t/* , de- 
note matrices on the link (n, e;)). The Us and Vs satisy 
(K) N = m N = 1, Vjpi = wU^V* (w = e 3 ^), i.e., 
Weyl's group relations, and matrices on different links 
commute. Zjv GFTs have been known for many years 
to be self-dual for N = 2,3,4, and it was conjectured 
that they are no longer self-dual for N > 5 14j. We can 
prove that these theories remain self-dual for all N, as 
the mapping of bonds 



V 2 

v n 

v 3 



AO 2 



A6l 



n—e 1 +e2 ' 



AO: 



n— ei+e3 > 



Ml 

AOl 



i t 



n — ei+e2+e3 
2 t 



V _L_ 

' n+e 3 
* n-\-e2 



(4) 



shows. Ujj is a new discrete symmetry of this problem, 
but Up = 1 up to a lattice translation. For large N, these 
gauge theories are known to display three phases, two 
of them connected through a confinement-deconfincmcnt 
phase transition [l7j]. The self-duality fixes the self-dual 
coupling g* at 4g* 2 K* = 1 [r|, which gives the exact self- 
dual coupling for every N (so far only known analytically 



for N = 2,3,4). On the other hand, it is shown in [15 1 
(using our approach) that the isotropic d+l = l + liV- 
state vector Potts model has a self-dual point at coupling 
J* given by precisely an equivalent relation 2K* = J*. 



Thus our results explain the puzzling fact 14J that the 
isotropic classical d + 1 = 1 + 1 TV-state vector Potts 
model and the d + 1 = 3 + 1 Z^v GFT share identi- 
cal self-dual relation: first, both bond algebras (though 
non-isomorphic) are based on the Weyl algebra, and ad- 
mit self-duality mappings; and second, both models have 
quantum couplings satisfying the equation in (l6| . The 
compactness of degrees of freedom (i.e., angular vari- 
ables), is required for a phase transitions to occur. On 
one hand, Polyakov [3] showed that compact QED dis- 
plays no phase transitions in 2 + 1 dimensions. On the 
other, we can show that in the limits TV — > oo, a — ► 0, Q 
reduces to the well know self-duality of vacuum QED in 
3+1 dimensions E i— ► B, B i— ► —75, which has no phase 
transitions. We argue that since the self-duality emerges 
only in 3 + 1 dimensions, it is important in triggering the 
phase transitions of these GFTs. So, the presence of both 
compactness and self-duality are crucial for the existence 
of a confinement-deconfincmcnt phase transition. 

In summary, we developed a unifying and systematic 
framework for dualities, providing a new perspective to 



unveil them: (sclf-)dualitics (exact or emergent) can be 
investigated as homomorphisms of bond algebras. The 
power of this algebraic approach was exploited to obtain 
new self-dualities of confining Abclian GFTs in 3 + 1 di- 
mensions, a new discrete symmetry of these theories, and 
their self-dual couplings analytically. We prove that the 
puzzling connections between these GFTs and some con- 
fining theories in 1 + 1 dimensions (vector Potts model) 
result from these two models having similar algebraic 
structures and self-dualities. Self-dualities are more eas- 
ily discovered as automorphisms of bond algebras (quan- 
tum) than as relations between partition functions (clas- 
sical). Furthermore, they can generate otherwise hidden 
symmetries. Known classical dualities derived in the lit- 
erature by Fourier transformation [io[ can be obtained by 
our technique. Thus this work hints at a deep connec- 
tion between operator algebra homomorphisms and the 
Fourier transform to be at the root of the equivalence 
between classical and quantum dualities. Our approach 
to (self-)dualities is applicable to any system, and clears 
the way for the development of approximation schemes 
that preserve these peculiar symmetries. 

E.C. thanks V. Lunts for helpful discussions. 
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Unified approach to Quantum and Classical Du- 
alities: Supplementary Material 

In this section, we present further applications of our 
method. Our approach provides a new perspective on all 
dualities in physics. With few exceptions, the study of 
(self-)dualities in classical and quantum systems has been 
dominated [l[ by techniques that amounts to a change 
of (classical) variables in partition functions. Our novel 
operator technique is different enough to require careful 
examination through several examples. 

Dualities of the extended Toric Code (TC) Model. The 
TC model Q was recently extended to include an exter- 
nal magnetic field h Q 

Hejc = -J X J2 A ° - J ^1 B p + ( 5 ) 

s p (ij) 

In Eq. ([5]), on each link of a square lattice there is an 
spin-1/2 operator, A s = Yij o~j S represents a product over 
the ^-component of the four spins that have s as a com- 
mon vertex, and B p ~ . o~j P is the product of the four 
cr| p operators that belong to the plaquette (□) p. The 
quantumphase diagram of -Hetc has been studied quite 
recently p|, 0] • It shows reflection symmetry relative to 
the line hx = h z , and a multi-critical point on this line 
as well 0, HI ■ The special role of the condition h x = h z 
can be understood in terms of a quantum self-duality. 
Drawing straight lines through the centers of the bonds, 
Hetc can we re-written in terms of spins lying on the 
vertices of another square at a 45 degree angle with the 
original one. On this lattice, the mapping 3>(er5) = cr5 + - 2 , 
$(cr5) = a~ + g 2 extends to a self-duality automorphism of 
Ah etc that exchanges J x with J z and simultaneously h x 
and h z . The reflection symmetry J x <-> J z and h x <-> h z 
relates to the Wegner duality and the more general 
self-duality of an Ising matter coupled gauge theory. On 
a cubic lattice, the action for the latter reads 

S = - J ViVijVj - K D ^ (6) 

{ij) □ 

where matter fields {ai} live at lattice sites i, while gauge 
fields rjij live on the links connecting sites i and j. This 
action is self-dual under 

J —> J = K, K -> K = — lntanh J, (7) 

which can be derived as a consequence of the self-duality 
of the ETC Model. This follows from the Euclidean 
representation of Eq. (|5|) given by the 3 (or 2+1)- 
dimensional action of Eq. ((6]) Q , with 

J = h z Ar, K = - lntanh^ At), (8) 

where At is the discretization step in the imaginary time 
direction. Inserting the duality crfj «-> afj into Eq. ©, 
we derive the gauge theory dualities of Eq. (J7J). The 



Wegner duality between the Ising model (K = 0) and 
the gauge theory (J = 0) is a particular case of Eq. (8]). 
The Ising matter coupled gauge system offers another ex- 
ample of an emergent duality. Take the system of Eq. ([6]) 
in the limit K — ► 00, that enforces a projection onto a 
space in which rjij = WiWj (uii = ±1). Setting Sj = aiWi, 
we obtain a general (i-dimensional rendition of the mat- 
ter coupled gauge theory of Eq. ((6]), the d-dimensional 
Ising model. The KW self-duality of the classical d = 2 
Ising model ( J — ► J = tanh _1 (exp[— 2 J])) appears as an 
emergent duality in the K — > 00 limit of the system of 
Eq. ([6]). In this limit, the bonds Zij = criTjijaj [8j satisfy 
the constraints Crj: Ylijen z ij = 1 f° r au plaquettcs □. 
In the projected subspace in which the constraints {Co} 
are satisfied, the algebras of the matter coupled gauge 
theory of Eq. © and the Ising model are identical. For 
finite K, the Ising matter coupled gauge system of Eq. 
(EI) 01 : is dual in d = 2 to an Ising model in a uniform 
magnetic field Q which does not obey the KW relations. 

Next, we introduce a model that is dual to the ETC 
model for arbitrary couplings. The dual Hamiltonian is 

H = -J^Urf- J z ^l4 

~h x J2 /''/'!<; • X! /'' /'}'/'.• ■ ( 9 ) 

{ij) 

In the particular case J z = h z = 0, this is the clas- 
sical Ising matter coupled lattice gauge theory. For a 
square lattice of N sites with periodic boundary condi- 
tions, there are matter fields of the type {fxf } and (2N) 
gauge fields For a dual system on a square lattice 

of N' sites, there are (27V') spin fields in the Hamiltonian 
of Eq. ([5]) . As the bond algebras in the two systems de- 
fined by the Hamiltonians of Eqs. (O, and ([9]) are the 
same then, ignoring additional global constraints stem- 
ming from boundary conditions, in computing the par- 
tition function Z, there is a one-to-one correspondence 
between the terms for the two dual models. If the only 
element with a non-vanishing trace is the identity oper- 
ator then, when N' = (3N)/2 then each of the terms in 
the expansion of Z will give an identical contribution. 
Finally, if h x = 0, a model dual to Eq. ([5]) is given by 

H = -j x ^ti-j,Y i uii*-h,Y l fiiij- 

i 1 (ij) 

The Hamiltonian of Eq. (fT0|) is determined by the N' 
variables {fxf }. This duality is exact if we scale the num- 
ber of variables accordingly: N' = 2N. 

Beyond S=|. The Blume-Emery-Griffiths (BEG) 
model Q -in its square lattice version (coordination 
z=A)- is given by an S=l Ising-like Hamiltonian 

PH BEG = -J2(JS?S* + K(S*S?) 2 ) + D]T(Sf) 2 .(ll) 



The BEG model was developed in the context of liquid 
3 He- 4 He mixtures. At particular values for the param- 
eters K = 3J, and D = 2zJ, this model becomes the 
isotropic q = 3 Potts model, f3H P = -2J£ (y) 6(Sf, S?), 
(S z = ±1,0). 

It can be shown, using the transfer matrix technique, 
that the d=l quantum BEG model is equivalent to 



V3 



(12) 



where the bonds are 



Viv 
2 



1 + V2Sf 



K 



D 



= 2 + s!si +x +-{s*s* +1 y-—{{sty + (s*y) 



V3 



J 



4 J 




with dual coupling constants B = 



2+A: 



In Xo 



2+A: 



In A 



3- 



D 



■HI 

AS 



In Ai 
2 



2+Af 



In A 2 



2+A; 



lnA 3 , and D ± = 



In A 2 



1-A2 

2+Af 



2+A 2 1 2+A 2 In A 3 . The parameters Ai,2,3 are given 

by Ai = 2e( K -T) s inh J, A 2 + A 3 = 2e^ K - ^ cosh J + 
1, and A 2 A 3 = 2(e( A '-*)coshJ - e~%). Finally, A Q = 
eS(l-Aa). 

The classical BEG system has a tricritical point, which 
is mapped to a quantum critical point of //qbeg- At 
K = 3J, D = 8 J (the Potts limit flO|). the bonds Ui and 
Vi satisfy the Hcckc relations iOiVi+\uji — u>i = Vi+iLUtVi+i — 

V l+1 , LJiUi — V%UJ l , V i U i V l - Vi = LUtViLOi - CJj , V t Vi = 

V3vi. Therefore, a quantum self-duality defined on gen- 
erators by Vi uJi, LOi — > Vi, fixes the self-dual point 
at 3B = 2J. Thus by algebraic means alone, we can 
determine the tricritical temperature T c from the dual 



coupling B: ksT c 



2. j 



in agreement with the re- 



ln(l+V3) ' 

suit obtained through other methods . 

Self- duality of vacuum Quantum Electrodynamics 
(QED). The classical (self-) duality of electromagnetic 
fields in the absence of sources (in vacuo) E — > B, B — * 
— E is one of the oldest examples of a self-duality avail- 
able. The bond algebra formalism gives insight into how 
this self-duality extends to non-compact vacuum QED. 
To see the connection with the self-duality we found 
for Zn lattice gauge theories, we quantize electromag- 
netism in the axial gauge, = 0. In this gauge, the 
classical fields in terms of the vector potential A arc 
E = — B — V x A. If we next apply canonical 

quantization to electromagnetism written in this gauge, 
we get vacuum QED in the form 



H 



d 3 x i(n 2 + (V x A) 2 ), 



together with 
[A m (x,t),Ili{y,t)] 
the canonical substitution ^ 



canonical commutation relations 
= i5 m i5(x — y), where we made 
-E — > II. Furthermore, 
from gauge invariance, we have that the subspace of 



FIG 

A 1 



1: The plaquette variable A6„ = A^ +ei — A% — A\ +E2 + 



physical states (gauge invariant states) of the full Hilbcrt 
space is specified by the so called Gauss constraint: 
V • II I physical) = (physically, this means we only keep 
those states on which V • E = 0). To investigate the 
bond algebra of this Hamiltonian, it is convenient to 
discretize the theory and consider it on a d = 3 cubic 
lattice of lattice spacing a. The Hamiltonian then reads 
(n = (n 1 , n 2 , n 3 )) 



^ 2 

=1 



where [A^,Il£j = i5 lk 5 n , m , and 

Ml = {Al +e2 -Al-Al +e3 +Al)/a 
M n = (A n+e3 — A n — A n+ei + A^)/ a 
= (A n+ei — A 2 n — A n+e2 + A\)/a 



are the discretized form of the components of V x A. We 
call the interaction terms A6 plaquette interactions, bor- 
rowing the terminology we used with Zjv gauge theories 
(see Fig. 1). The first thing to notice is that only four 
plaquette variables AO have non-trivial commutation re- 
lations with any given momentum component 11^ at site 
n and similarly, only four II fields have non-trivial com- 
mutation relations with any given A9 l n . Thus we have a 
good case for an automorphism that exchanges II «-> Ad. 
In fact, one such automorphism is given by 



K 
n 3 h 



>n r \ 
>n 2 
ri 3 : 



A6l ^ AB 2 n 
ABl » Adl 



n 1 

n— ei+e 2 +e3 



= A0 n _ ei+e2 , 

LAO n-e 1 +e 3 > 



-K+e 3 (13) 

n 3 

11 n+e 2 ■ 



The geometry of this mapping is clarified in the Figs. 2 
and 3. A moment's reflection makes it clear that this 



7 



C2 




e3 



FIG. 2: Schematic of the bond algebra mapping (|13[1 . Associ- 
ated with each lattice site n, there are three II fields and three 
A9 plaquettes. The A9 plaquettes at site n map by (| 1 3fl to 
displaced lis, each colored plaquette mapping to the —II at 
the correspondingly colored site. 

mapping is nothing other than E i— > B, B — E, the 
quantum descendant of the classical electromagnetic self- 
duality. As explained in the main body of our paper, the 
self-duality mapping serves the double purpose of estab- 
lishing the existence of a self-duality and defining the 
dual variables. We have explicitly shown the dual vari- 
ables above, naming them II and A. The dual vector 
potential is defined implicitly by the above relations, as, 
for instance, A0* = (A 3 n+e2 - A 3 n - A 2 n+ea + A 2 n )/a. 

It is interesting to compute explicitly the dual variables 
in the continuum limit a — > 0, paralleling the discussion 
of dual variables for the scalar field given in the main 
body of this paper. In this limit, Q13p implies the rela- 
tions 

n(a;,i) = V x A(x,t) 
V x A(x,t) = -Tl(x,t) 

between the initial and the dual operator variables. Thus 
II is already explicitly given in terms of A. We need to 
solve the second relation for A. It is not difficult to check 



that, on physical states (on which V ■ II vanishes), 
A(x,t) = -— V x / d 3 w , , . 

V ' 47T J \X~W\ 
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FIG. 3: The effect of the exchange duality of (|13[1 on the three 
II fields at site n. 



